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Dynamics of ordering in Ising model, following quench to zero temperature, have been studied 
via Glauber spin-flip Monte Carlo simulations in space dimensions d = 2 and 3. One of the primary 
objectives has been to understand phenomena associated with the persistent spins, viz., time decay 
in the number of unaffected spins, growth of the corresponding pattern and its fractal dimensionality, 
for varying correlation length in the initial configurations, prepared at different temperatures, at and 
above the critical value. It is observed that the fractal dimensionality and the exponent describing 
the power-law decay of persistence probability are strongly dependent upon the relative values of 
nonequilibrium domain size and the initial equilibrium correlation length. Via appropriate scaling 
analyses, these quantities have been estimated for quenches from infinite and critical temperatures. 

The above mentioned dependence is observed to be less pronounced in higher dimension. In addition 
to these findings for the local persistence, we present results for the global persistence as well. 
Further, important observations on the standard domain growth problem are reported. For the 
latter, a controversy in d = 3, related to the value of the exponent for the power-law growth of the 
average domain size with time, has been resolved. 


I INTRODUCTION 

Kinetics of phase transitions [l|-[J remains an ac¬ 
tive area of research for several decades. In this area, 
typically one is interested in the nonequilibrium dy¬ 
namics related to the evolution of a system to a new 
equilibrium state, having been quenched from a con¬ 
figuration prepared outside the coexistence curve to 
inside it, via the variation of temperature (T), pres¬ 
sure, etc. In this work, our focus is on the paramag¬ 
netic to ferromagnetic transition When a system 
is quenched, via variation of T, from the paramag¬ 
netic phase to ferromagnetic on^ domains rich in like 
spins form and grow with time [2[ . Aspects that drew 
attention of researchers, in this problem, are under¬ 
standing of domain patterns [H, growth of domains 
[H, aging properties of the evolution [gI-Q, as well 
as the pattern (and corresponding dynamics) exhib¬ 
ited by atomic magnets (or spins) that did not change 
orientation till time t, referred to as persistent spins 
iihll- This work deals with issues related to do¬ 
main growth and persistence. 

During the process of ferromagnetic ordering (where 
the order parameter is a nonconserved quantity), the 
average domain size, i, increases as [H 

ir^t^, (I) 

where a, the growth exponent, may have dependence 
upon system dimensionality (d) based on the order- 
parameter symmetry. This growth occurs via motion 
and annihilation of defects, facilitated by change in 
orientation of the spins, Si, the subscript i being an 
index related to an atom or spin, typically considered 
to be located on a regular lattice. In this work we 
study the spin-I/2 Ising model, to be defined later, 


for which defects are the domain boundaries. In this 
case, Si is a scalar quantity which gets affected only 
via (complete) flipping or change in sign. For this 
model, the theoretical expectation for a is same in 
both d = 2 and 3. 

The persistence probability, P, defined as the frac¬ 
tion of unaffected spins, typically decays as Q 

P - (2) 


where 0 is expected to have dependence upon d. 
The persistent spins exhibit interesting fractal pattern 
with dimensionality [l7| df whose dependence upon 
0 will be introduced later. Unless mentioned other¬ 
wise, all our results on this issue correspond to lo¬ 
cal persistence, probability for which, as already men¬ 
tioned, is calculated by counting unaffected “micro¬ 
scopic” spins. There has also been interest in the cal¬ 
culation of such probability by dividing the system 
into blocks of linear dimension and counting the 
persistence of coarse-grained or block spin variables 
|l4 llbf . In the limit -£5 ^ a, the microscopic lat¬ 
tice constant, such block persistence probability, P5, 
will correspond to P, the “local or site persistence” 
probability. On the other hand, for ^5 ^ 00, one 
obtains “global persistence” probability, further dis¬ 
cussion and results for which will be presented later. 

For Ising model, values li [1^ . [Toj of a, 0 and df 
are accurately estimated via Monte Carlo (MC) sim¬ 
ulations, in d = 2, for quenches from initial temper¬ 
ature Ti = 00 to the final value 77 = 0. It is rea¬ 
sonably well established [H, 0, [2^ that, in this 

case, the values of a, 0 and df are 1/2, 0.225, and 
1.58. However, the conclusions, if exists, on the corre¬ 
sponding numbers for d = 3 are questionable iGl-lisj . 
Recent focus, on the other hand, for persistence as 
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FIG. 1. Log-log plots of local persistence probability, P(t), 
vs t, for quenches (of the Glauber Ising model) from differ¬ 
ent values of initial temperature Ti {>Tc, the critical tem¬ 
perature), to the final value T/ = 0. All results correspond 
to space dimension d = 2 and square lattice, with linear 
dimension of the square box being L = 2048, in units of 
the lattice constant a. The lines represent various power- 
law decays, values of the exponents being mentioned in 
appropriate places. 

well as for other aspects of coarsening, has been on 
[idl-liBI. [29l - [^ quenches from temperatures providing 
large equilibrium correlation length In this context, 
in a recent work [2^ , we have explored the initial cor¬ 
relation dependence of a and 0. Our observation was, 
while a is insensitive to the variation of Ti (at least in 
d = 2), P (and thus 0) is strongly influenced by the 
choice of the latter, viz., we obtained for d = 2 and 
3, d = dc — 0.035 and 0.105 for Ti = Tc, the criti¬ 
cal temperature (see Fig. [T]for d = 2). The numbers 
quoted above are significantly different from those for 
Ti = 00 . For intermediate temperatures, as seen in 
Fig. m two step decays can be noticed. A slower 
decay was observed for £ < ^. The corresponding ex¬ 
ponent 0i approaches dc as ^ ^ oc, i.e, when Ti Tc. 
For ^ ^ C, behavior consistent with Ti = 00 was ob¬ 
tained. This implies, dynamics of the spins is strongly 
influenced by the relative values of nonequilibrium do¬ 
main length £ and the equilibrium correlation length 
^ in the initial configuration. The overall time decay 
of P, for all Ti, was empirically constructed to be [2J 


(d — 0i)/a^ 2 , whereas Co and Ci are dimension 

dependent constants. 

An extension of a study 0 (via a different model 
in d = 1 ) predicts 

df = d — zO^ (5) 

where 2 ;, to be more formally defined later, is a dynam¬ 
ical exponent related to the growth of the persistence 
pattern. From previous studies even though it 

has been reported that the decay of P is disconnected 
with the growth of £, z and a may be related. Nev¬ 
ertheless, since such a connection is unclear, to gain 
knowledge about the variation of d/, as a function of 
Tj, estimation of z is needed. Even if such a connec¬ 
tion exists, as mentioned, the value of a in d = 3 is 
not unambiguous. In this dimension, the theoretically 
m expected value of 1 / 2 ) disagrees with some 

computer simulations [ 1 ^ which report numbers close 
to 1/3. This difference can possibly [ 2 ^ be due to 
long transient period. Thus, lengthy simulation runs 
with large systems are needed. It will be interesting to 
see if such long simulation, luck favoring, can provide 
the theoretically expected value. If yes, in that time 
regime, do we see change in other quantities as well? 

In this work, our objective thus, is to estimate d/, 
2 ;, a and d, for T^ = 00 and Ti = Tc, in space dimen¬ 
sions 2 and 3, for quenches to Tf = 0. For the ease 
of reading, in TABLE 1 we provide a list of values of 
these quantities, obtained from computer simulations. 
While the ones with asterisks, to the best of our knowl¬ 
edge, will be calculated (or the simulation results will 
be shown to be consistent with those theoretical ex¬ 
pectations) for the first time, the numbers appearing 
inside the parentheses are improvements over the ex¬ 
isting ones that appear outside. We will start present¬ 
ing results with the objective of calculating df. Other 
quantities will be needed for this purpose and will be 
estimated in due course. 


TABLE 1. List of some nonequilibrium exponents for Ising 
model. 


Gase 

a 

2: 

0 

df 

d = 2 , Ti = 00 

1/2 

2 

0.225 

1.58 (1.53) 

d = 2 ,Ti = Tc 

1/2 

2* 

0.035 

1.92* 

d = S, Ti = 00 

1/3 (1/2) 

2* 

0.18 (0.15) 

2.65* 

d = 3 ,Ti = Tc 

1/2 

2* 

0.105 

2.77* 


= A 


g(x) + : 




(3) 


with 


9(x) = 


Co 

1 + Cix^ ’ 


(4) 


where A is the amplitude of the long time decay, (j) = 


In Eig. [2] we show persistence snapshots for Ti = oc 
and Tc, both from t = 10 ^ Monte Carlo steps (MCS), 
this time unit to be defined soon, for d = 2 Ising 
model. It is clear that the patterns are different and 
so, different values of df are expected. In Eig. O 
snapshots from an intermediate temperature T^ = 2.4 
(> Tc), for d = 2, are presented. The first frame corre¬ 
sponds to a time falling in the slower decay regime of 
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Ti = oo Ti = Tc 



FIG. 2. Snapshots of the persistent spins are shown for 
quenches from = oo and Tc, to Tf = 0. The results 
correspond to d = 2, L = 2048 and t = 10^ MGS. In 
both the cases only parts of the boxes are shown where 
the persistent spins are marked in black. 


FIG. 3. Snapshots of the persistent spins from different 
times, mentioned on the figure, are shown for Ti = 2.4 
and T/ = 0. Other details are same as Fig. [21 


where J is the interaction strength (> 0) and < ij > 
implies interaction among nearest neighbors. The val¬ 
ues of Tc for this model in d = 2 and 3 are respectively 
[ 3 ^ 2.21J/kB and 4.51 ks being the Boltz¬ 

mann constant. 

Kinetics in this model was introduced via Glauber 
spin-flip mechanism [sH, [33[ . In this MC approach, a 
trial move consists of changing the sign of a randomly 
chosen spin. Since our quenches were done to Tj = 0, 
a move was accepted only if it had reduced the energy. 
Needless to say, initial configurations were prepared at 
nonzero T values. In that case, the Metropolis crite¬ 
rion for the acceptance of a move was implemented 
via appropriate calculation of the Boltzmann factor 
[in and its comparison with a random number, rang¬ 
ing between 0 and 1, whenever the move brought an 
increment in the energy. For preparation of initial 
configurations at temperatures very close to Tc, in 
addition to the Glauber mechanism, we have applied 
Wolff algorithm [13 as well, which facilitates faster 
equilibration. Time, in our simulations are measured 
in units of MCS, each MCS consisting of steps, L 
being the linear dimension of a square or cubic box. 
Periodic boundary conditions were applied in all di¬ 
rections. Final results are presented after averaging 
over multiple initial realizations, the number ranging 
from 20 to 70. In d = 2 all results are for L = 2048. 
In d = 3, the results for Ti = 00 are for L = 512 and 
for Ti = Tc^ we presented results from L = 400 and 
256. 




Fig. H] (for the corresponding temperature), whereas 
the second one is from the faster decay regime, imply¬ 
ing i ^ The earlier time snapshot resembles the 
Ti = Tc picture of Fig. [2] and the second one has sim¬ 
ilarity with Ti = 00 pattern. This justifies our focus 
only on these two limiting initial temperatures with 
^ = 0 and 00 , rather than exploring a wide tempera¬ 
ture range, to accurately quantify df and z. 

The rest of the paper is organized as follows. In the 
next section we describe the model and method. Sec¬ 
tion III provides a brief overview of an earlier work. 
Results are presented in section IV. Finally, section V 
concludes the paper with a brief summary and out¬ 
look. 


Ill AN OVERVIEW OF THE BACKGROUND 
ON FRACTALITY OF PERSISTENCE 
PATTERN 


In this section we provide a discussion on the the¬ 
oretical background for fractality of the structures 
formed by persistent spins, following the work by 
Manoj and Ray 

From a density correlation function, T(r, t), 
isotropic in an unbiased system, total mass or number 
of particles in a circular or spherical (depending upon 
dimensionality) region of radius R can be obtained as 


M{RR) 


.R 

/ D{r,t)i 
Jo 


.d-1 


dr. 


( 7 ) 


II MODEL AND METHOD 

As already mentioned, we study the Ising model Q, 
on square or simple cubic lattice systems, depending 
upon the dimensionality, with nearest neighbor inter¬ 
actions. The Hamiltonian for the model is given by 


r (= |r|) being the scalar distance of a point in that re¬ 
gion from the central one. An appropriate correlation 
function in the present context is 


T(r,t) 


(/)(ro,t)p(ro +f,t)) 

{piro,t)) 


(8) 


H = -JJ2 SiSj; Si = ± 1 , 

<ij> 


with p being unity at a space point if the spin there 
did not flip till time t and zero otherwise. The average 
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order parameter for the persistent pattern is 

{p{r,t)) = ^ = P{t). (9) 

This being a nonconserved (time dependent) quan¬ 
tity and, since, in the definition of D{r,t), the av¬ 
erage value is not subtracted from p, decorrelation 
here means, decay of D{r^t) to a “non-zero” value 
{=P(t))^ for t < cxo. The distance, ip{t), at which 
D{r, t) reaches this plateau is the characteristic length 
scale of the pattern. In that case, there may exist scal¬ 
ing of the form 

^ fir/1,). ( 10 ) 

For X {= r/Ip) > 1, / should be unity. On the other 
hand, for fractal dimension df and x < 1, one should 
have 

fix)^x‘^f-'^, (11) 

since 

Mr^X^P (12) 

Considering that P(t), the plateau value, decays in 
a power-law fashion, a power-law behavior of f{x) is 
indeed expected, once scaling is achieved. A continu¬ 
ity, at r = in such a situation demands 

l{df-d)lz ^ 

providing Eq. m, where 2 ) is the dynamic exponent 
characterizing the growth of the persistence pattern, 
mentioned before, as 

ip ~ i'/". (14) 

For this model, as mentioned, value of a has been 
estimated 0 for various Ti values in d = 2. However, 
a priori it is unclear whether there is a general validity 
of the relation 

za = l. (15) 

Then it is necessary to calculate both z and d, for 
correlated and uncorrelated initial configurations, to 
validate Eq. ([5j). On the other hand, as already men¬ 
tioned, the value of a is ambiguous in d = 3. 

IV RESULTS 

In Fig. m we show D{r,t) as a function of r for 
Ti = 2.4, from two different times, mentioned on the 
figure, for d = 2. As expected, the correlation function 
decays to different constant value, P(t), at different 
length £p^ for different times. Before decaying to the 



FIG. 4. Density correlation functions, D(r, t), related to 
the persistent spins, are plotted vs r. Results are presented 
from two different times, for Ti — 2.4 and T/ = 0. The 
system dimensionality is d = 2 and value of L is 2048. 


plateau, the early time data appear to obey a power- 
law. The later time data, for smaller r, follows the 
same power-law before crossing over to another, faster, 
power-law decay. This implies, there exist two length 
scales in the problem, below and beyond the equilib¬ 
rium scale Inside the larger structure, the small 
length scale structure remains hidden, which will be¬ 
come irrelevant in the long time limit. For ^ = oc, i.e., 
Ti = Tc, however, the latter will be the only structure 
and remain for ever. The exponent for large r, for 
Tc < Ti < oo and t ^ 0, should be related to the df 
value for Ti = oo case whereas, in case of small r, the 
exponent should be connected to df for Ti = Tc case. 
Below we focus on these two cases, i.e., Ti = oo and 
Ti = Tc^ separately, first for d = 2, followed by d = 3. 
As need occurs, we will present results related to a, 0 
and 2 ). 

In the main frame of Fig. [5] we present a scaling 
exercise lillii for D{r) where we have plotted f{x) 
as a function of x, using data from different times af¬ 
ter quench, for Ti = Tc and d = 2. Scaling appears 
good and gets better with the progress of time. On 
this log-log plot, look of the data appear, before de¬ 
caying to unity, linear, implying a power-law decay. 
The exponent appears to be 0.09. In the inset of 
this figure, we show analogous exercise for Ti = oo. 
Even though this case in this dimension was studied 
by Jain and Flynn , for the sake of comparison and 
completeness, we present it here from our own simu¬ 
lations. In this case, the exponent for the power-law 
decay appears consistent with 0.45. Then, in d = 2, 
for Ti = oc, the fractal dimensionality is 1.55 and for 
Ti =Tc, the number is 1.91, if Eq. dU is valid. 

In Fig. [6] we show the plots of vs t in d = 2, 
for (a) Ti = oo and (b) Ti = T^ on log-log scales. In 
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FIG. 5. Scaling analysis of D{r, t) for the d = 2 Ising 
model, with Ti = Tc and T/ = 0, where f{x) is plotted vs 
X = r/ip, using data from different times after the quench, 
on log-log scale. The solid line corresponds to a power-law 
decay with an exponent 0.09. Inset: Same as the main 
frame but for Ti — oo. The solid line here has the power- 
law decay exponent 0.45. The value of L is 2048 for all 
the results. 


both the cases the data appear consistent with 2 : = 
2, validating Eq. (note that a is established to 

be 1/2 in d = 2). Nevertheless, we intend to make 
more accurate quantification. For this purpose, in the 
insets of these figures we have shown instantaneous 
exponents (dash-dotted lines), calculated as [35| 


1 d\n£p 
Zi dint ^ 


(16) 


YS 1/ip. In both the cases we obtain the value of 
z via linear extrapolation (see the consistency of the 
simulation data with the solid line) to = 00 . For 
= 00 , from this exercise, we quantify 2 : = 2.15 and 
for Ti = Tc, we obtain the number 2 : = 2.02 (see the 
dashed horizontal lines). These numbers, in addition 
to verifying Eq. dm, are also consistent with the 
numbers obtained via least square fitting of the ip vs 
t data to the form 


= (17) 


where i\ and Az are positive constants. This consis¬ 
tency may imply, early time corrections to the expo¬ 
nents are insignificant. Note here that, in absence of 
any correction, one expects 


1 


i^ 

1 - ^ 


(18) 


a linear behavior of 1 / 2 :^, when plotted vs 1/ip, with 
slope —i^/z. A positive slope in both the insets is 



FIG. 6. (a) Log-log plot of persistence length scale, ip^ as a 
function of t, for d = 2 Ising model, following quench from 
Ti = 00 to T/ = 0, with L — 2048. Inset: Instantaneous 
exponent, Zi, obtained using the data in main frame, is 
plotted vs. l/ip. (b) Same as (a), but here Ti = Tc. The 
solid lines in the main frame of both (a) and (b) correspond 
to power-law growths with exponent 1/2. In the insets, the 
horizontal dashed lines correspond to our estimates for z, 
whereas solid lines are guides to the eyes. 


due to the fact that we have presented inverse of the 
quantity discussed in Eq. m- Using these values of 
2 ), and numbers for 0, mentioned earlier, in Eq. 
we obtain df = 1.93 for T^ = Tc and df = 1.51 for 
Ti = 00 . These values, within computational errors, 
are consistent with the conclusions from Fig. [5l Next 
we present results from d = 3. 

In d = 3, we start by presenting results for the 
growth of i, considering the controversy 14, 23, 2^ on 
the value of a discussed above. In d = 2, we avoided 
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presenting results on this aspect with the understand¬ 
ing that the issue there is well settled. Nevertheless, 
in the context of global persistent decay, we will make 
indirect conclusion about it. Here note that the esti¬ 
mation of ^ was done from the first moment of domain 
size distribution, p{£d,t)^ as 


e{t) = j edp{£d,t)d£d, (19) 

where id is the distance between two domain bound¬ 
aries in a particular direction. The main frame of 
Fig. [71^a) shows a plot of i vs t, on log-log scale, for 
quenches of the d = 3 Ising system from = oc to 
Tf = 0. There exists an intermediate time regime, 
extending over more than two decades, during which 
the simulation data show consistency with an expo¬ 
nent a = 1/3, in agreement with previous results 
M- However, as discussed and a trend demonstrated 
in Ref. [1^, the discrepancy in the earlier reports 
from the theoretical number 1/2 can be due to long 
transient. Thus, long simulation runs with large sys¬ 
tems are necessary. We have simulated a system with 
L = 512 for time longer than any of the previous 
works, to the best of our knowledge. Indeed, it ap¬ 
pears that the long time behavior, over the longest 
time decade in the presented time range, is consistent 
with a = 1/2. In the inset of this figure we show the 
instantaneous exponent 


(Xi — 


dlni 

d\nt^ 


( 20 ) 


as a function of 1/i. This provides an accurate pic¬ 
ture, the long time exponent being within 5% of the 
theoretical value. One may then ask, is the value of 0 
going to change, beyond this crossover time? Even if 
0 has no dependence on the value of a, such a change 
may still occur. Note that conclusion on the value 
of in earlier works 3, 241, were drawn from runs 
shorter than this. Indeed, a jump in 6^^, calculated 
from 


o^ = - 


d\nP 
dint ’ 


( 21 ) 


occurs (see Fig. [Tfb) and the corresponding inset) 
from an early time value of 0.18 to 0.15. This 
may, of course, be due to statistical or other reasons. 
However, since the jump in ai occurs around the same 
time as the one for Oi and fluctuation is seen around 
stable mean values, in both and we accept this 
as the correct number for 0 in the asymptotic time 
limit. 

Whether due to lattice anisotropy Q oi* anything 
else, the solution to overcome such long transient is 
certainly related to being able to access large length 
scales. For Ti close to Tc, since this is automatically 
the case, due to large initial correlation, we expect an 


enhanced value of a from early time. Corresponding 

vs t data are presented in Fig. [51 On the log-log 
scale, this data set shows consistency with a = 0.45. 
Here we mention that study for Ti = Tc has addi¬ 
tional problems related to longer equilibration time 
at the initial temperature and stronger finite-size ef¬ 
fects during the nonequilibrium evolution [ 2 ^. The 
latter remark can be appreciated from the plot in Fig. 
[8] where a bending of the data set (from the power- 
law behavior) is visible from t = 10^. This should be 
compared with the corresponding data in Fig. [Tl^a) 
for Ti = oc. Thus, accessing very large length scales, 
without finite-size effects, for Ti = Tc, is extremely 
difficult. The P vs t data, shown in the inset of Fig. 
ISl exhibit consistency [2^ with Oc — 0.105. Since, a 
is very close to 1/2 already, we do not expect much 
change in Oc even in true asymptotic length or time 
limit. The calculations of ai and Oi in this case pro¬ 
vide numbers consistent with the ones quoted above. 

Next, we come back to the issue of fractality. For 
d = 3 Ising model, our results in this context are en¬ 
tirely new irrespective of the value of Ti. 

Fig. ini^a) is analogous to the inset of Fig. ^{Ti = 00 
results for f{x)) but for d = 3. The corresponding 
f{x) vs X scaling plot for = Tc in d = 3 is presented 
in Fig. mb). Again, for both Ti = 00 and Ti = T^ 
good data collapse are obtained for results from dif¬ 
ferent times, in these scaling plots. For both values 
of Ti, we have used data sets lying in the time ranges 
that provide consistency with the expected theoreti¬ 
cal number for a. In the relevant region, the Ti = Tc 
results have power-law decay with exponent 0.24. In 
case oi Ti = 00 , the value of this exponent is approxi¬ 
mately 0.38. These numbers imply df = 2.76 and 2.62 
for Ti = Tc and Ti = oc, respectively. 

In Fig. [To] we show ip vs t plots from d = 3 for (a) 
Ti = 00 and (b) Ti = Tc. In the long time limit, the 
results, in both (a) and (b), appear consistent with 
growth having z = 2. This is in agreement with Eq. 
(fT5]) . From the log-log plot for Ti = oc, like i vs t, 
a long time transient is clearly visible. To quantify z 
more accurately (in the t ^ 00 limit), for both Ti = oc 
and Ti = Tc, we have shown the instantaneous expo¬ 
nents, vs 1/ip, in the insets. From there, we extract 
z = 2.1 ioT Ti = oc and 2.15 for Ti = Tc. Along- 
with the above mentioned numbers for z, using the 
values of 0 for quenches from Ti = Tc and Ti = oc, we 
obtain df cz 2.78 and 2.68. These numbers are consis¬ 
tent with those obtained from the scaling plots in Fig. 
m providing higher confidence on our estimation of 0 
from long time limit, for Ti = oc. An interesting exer¬ 
cise here would have been to plot Ziai vs t. However, 
a constant value of unity cannot be obtained because 
of the fact that i and ip have different initial off-sets. 
This is evident from the pictures in the insets of Fig. 
[3(a) and Fig. [TOf ab While for the time dependence 
of i, a long transient with a cz 1/3 is visible, this is 
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FIG. 7. (a) Log-log plot of ^ vs t, in 0 ? = 3, for = 00 . The 
solid lines correspond to different power laws, exponents 
for which are mentioned. The inset shows instantaneous 
exponent as a function of Xjt. The dashed horizontal 
lines represent exponent values 0.36 and 0.48. (b) Log-log 
plot of P(t) vs t, for d = 3 and Ti — 00 . The inset shows 
instantaneous exponent Oi vs l/l. Horizontal dashed lines 
are for 0 = 0.176 and 0.15. 


not so for the time dependence of Thus, because of 
the off-set related reason mentioned above, Ziai = 1 
is expected to be valid only in the t ^ 00 limit. 


Finally, we turn our attention to the block persis¬ 
tence which was introduced by Cueille and Sire |14|. 
The corresponding probability P 5 , as already men¬ 
tioned, is related to the change in the order-parameter 
variable obtained by coarse-graining the site or micro¬ 
scopic spin variables over a block of linear size -£ 5 . It 
is expected that the decay of this probability will be 
significantly slower than the site or local persistence 



FIG. 8. Log-log plot of i vs t, for d = 3, L = 400 and 
Ti — Tc. The solid line corresponds to a power-law growth 
with exponent 0.45. The inset shows a plot of P vs t, on 
log-log scale. The solid line there represents a power-law 
decay with exponent 0.105. 


probability, to which the former should cross over only 
for £ > This two time-scale behavior is desirable 
by considering that, in the early time regime, a slower 
decay is forced by the fact that a sign change in block 
spin variable happens only when £ becomes compara¬ 
ble to £b and in the large £ limit, the blocks effectively 
appear as sites. It is expected then that a scaling 
should be obtained as 

( 22 ) 


where Oq is the exponent of the early part of the decay 
or global persistence exponent in the sense that when 

^ 00 , this is the only exponent. In d = 2, we will 
see that the best scalings, in accordance with Eq. (j22j), 
are obtained for a = 1 / 2 , irrespective of the value of 
Ti. In d = 3, on the other hand, due to long transient 
in the dynamics, we avoid presenting these results. 

In addition to the above mentioned understandings, 
calculation of persistence probability via such block¬ 
ing may have advantage for quenches to nonzero tem¬ 
perature. Note that for Tf ^ 0, thermal fluctuation 
from bulk of the domains affects the calculation when 
done via standard method. Considering that domain 
growth occurs essentially due to spin flips along the 
domain boundaries, in the calculation of P, dynam¬ 
ics inside the domains needs to be discarded. In a 
method, prescribed by Derrida [l 2 |, this is done by 
simulating an ordered system, alongside the coarsen¬ 
ing one, and subtracting the common flipped spins, 
identifiable as the bulk flips, between the two systems, 
from the total, thus sticking to the effects of only the 
boundary motion. In the block spin method, if is 






























FIG. 9. (a) Scaling function f{x) is plotted vs x, for d = 3, 
Ti = oo and T/ = 0, using data from few different times. 
The solid line has a power-law decay with exponent 0.38. 
(b) Same as (a) but, instead of = oo, we present data for 
Ti = Tc. Here the solid line has power-law decay exponent 
0.24. The results were obtained for simple cubic lattice 
with L = 256. 


significantly larger than ^ at Tj, thermal fluctuations 
will not alter the sign of block spins and in the large i 
{> £b) limit, as previously stated, one expects the de¬ 
cay to be consistent with local persistence. This saves 
computational time for simulating the additional sys¬ 
tems with ordered configurations. 

In Fig. nra) we show Pb vs t plots from d = 2, for 
a few different values of ib and Ti = Tc. It appears, as 
discussed, there exist two step decays and crossover to 
the faster (consistent with the local persistent decay) 
one is delayed with increasing ib. 

In Fig. [TTI bi we show a scaling exercise using the 



FIG. 10. (a) Log-log plot of persistence length scale, ip, 
as a function of t, for d = 3 Ising model, following quench 
from Ti = oo to T/ = 0, with L = 512. Inset: Instanta¬ 
neous exponent, Zi, obtained using the data in main frame, 
is plotted vs. 1/ip. (b) Same as (a), but here T = Tc and 
L = 256. The solid lines in the main frame of both (a) 
and (b) correspond to power-law growths with exponents 
mentioned there. The horizontal dashed lines in the insets 
correspond to our estimates for z and the solid lines there 
are guides to the eyes. 

data of Fig. m a) where we have plotted h{x) vs 
. For obtaining collapse of data, we have ad¬ 
justed do and a. The value of a used here is 0.49, that 
provides the best collapse. This number is certainly 
consistent with 1/2, within numerical error. Early 
time behavior corresponds to global persistence with 
do = 0.002 and the late time behavior is consistent 
with our previous estimation of dc — 0.035, for the site 
persistence probability. In the inset we have shown 
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FIG. 11. (a) Plots of block persistence probabilities, vs 
t, from different values of £b, for Ti = Tc, in d = 2. 
(b) Scaling plots of the persistence probabilities in (a). 
The scaling function h{x) is plotted, on a log-log scale, vs 
X = . In the inset of (b) we present similar scaling 

plot for Ti = oo. System sizes correspond to L = 2048. 
Various power-law decays are shown by solid lines with 
the exponent values being mentioned next to appropriate 
lines. 


corresponding scaling results for T^ = oc, for which 
Oq and 0 values (mentioned on the figure) are consis¬ 
tent with previous findings Q. The value of a that 
provides the best collapse here is 0.5. Note that in 
our earlier work such independence of a from Ti was 
directly (from the analysis of ^ vs t data) checked for 
this dimension. 


V CONCLUSION 

We have presented results for coarsening dynam¬ 
ics in Ising model, with nonconserved order param¬ 
eter, from space dimensions d = 2 and 3. The re¬ 
sults include domain growth law and persistence, for 
quenches with initial configurations of varying correla¬ 
tion length While presented results for persistence 
are mostly related to local order parameter ISlii , for 
the global case Ullli we have obtained new exponent 
for quench from initial temperature T^ = Tc, in d = 2. 
For local persistence, our results are summarized in 
the next paragraph. 

A central objective of this paper has been to iden¬ 
tify the differences in the patterns formed by persis¬ 
tent spins when systems are quenched from T^ = oo 
and Ti = Tc^ to the final temperature Tj = 0. For 
both the cases, corresponding fractal dimensionalities 
d/, as well as the exponent z, related to the growth 
of the persistent pattern, have been obtained in var¬ 
ious dimensions. A scaling law connecting d/, d, z 
and 6>, predicted by Manoj and Ray 0 , has been 
observed to be valid, irrespective of the values of d 
and Ti. Combining various methods, we quote, for 
Ti = oo, 

df = 1.53T0.02, d = 2, ^ ^ 

^ 23 

d/= 2.65T0.03, d = 3, ^ ^ 

and for Ti = Tc, 

df = 1.92T0.02, d = 2, ^ ^ 

df = 2.77 ±0.02, d = 3. ^ ^ 

On the standard domain growth problem, it is 
shown that the values of a in both dimensions for 
all initial temperatures are consistent with the theo¬ 
retical expectation a = 1/2. This number describes 
the growth of the persistent pattern as well, validating 
Eq. ([T5|) and confirming that domain growth occurs 
essentially due to dynamics of spins along the domain 
boundaries. This resolves a controversy in d = 3 for 
which some previous computer simulations reported 
a = 1/3. As mentioned in Ref. [26|, this discrepancy 
must have been due to lack of data for extended period 
of time. Long simulations in our work, in addition to 
resolving this controversy, corrects the value of 0 as 
well in this dimension. 

In future we will address similar issues for conserved 
order parameter dynamics, including aging phenom¬ 
ena. For both conserved and nonconserved dynamics, 
scaling properties and form of the two-point correla¬ 
tion function will be an important problem for the 
case of correlated initial configurations. 
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